Long Range Bond-Bond Correlations in Dense Polymer Solutions 
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The scaling of the bond-bond correlation function Pi(s) along linear polymer chains is investi- 
gated with respect to the curvilinear distance, s, along the flexible chain and the monomer density, 
p, via Monte Carlo and molecular dynamics simulations. Surprisingly, the correlations in dense 
three dimensional solutions are found to decay with a power law Pi{s) ~ s~" with uj = 3/2 and the 
exponential behavior commonly assumed is clearly ruled out for long chains. In semidilute solutions, 
the density dependent scaling of Pi{s) « (^/s)"" with ujo = 2 — 2^ = 0.824 (i/ — 0.588 being 
Flory's exponent) is set by the number of monomers g{p) contained in an excluded volume blob of 
size ^. Our computational findings compare well with simple scaling arguments and perturbation 
calculation. The power-law behavior is due to self-interactions of chains caused by the chain con- 
nectivity and the incompressibility of the melt. This study suggests a careful reexamination of the 
operational definitions used for the experimental determination of the persistence length. 
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In this Letter we study the correlations of the direc- 
tions of bonds alonff polymer chains in semidilute solu- 
tions and melts 0,0,y|- We focus on flexible monodis- 
perse chains of N monomers (cf. Fig. ^ under good 
solvent conditions in three dimensions {d = 3) where 
both the bond length I and the excluded volume screen- 
ing length ^ are always much smaller than the chain 
end-to-end distance Re- In principle, once the bond-bond 
correlations are computed all other conformational sin- 
gle chain properties can be derived. Importantly, being 
the (second) derivative of the spatial distances along the 
chains, they allow us to probe directly — without trivial 
ideal contributions — the non-gaussian corrections pro- 
posed recently [^. As we shall see, these corrections are 
crucial to make the description of dense polymer systems, 
first proposed by Flory |3| and later corroborated by Ed- 
wards fully self-consistent. 

The bond-bond correlation function Pi (s) is generally 
believed to decrease exponentially Q . This belief is based 
on the few simple single chain models which have been 
solved rigorously 3, 6] and on the assumption that all 
long range interactions are negligible on distances larger 
than ^ due to the screening mechanism described by Ed- 
wards 0, 0. Hence, only correlations along the back- 
bone of the chains are expected to matter and it is then 
straightforward to work out that an exponential cut-off 
is inevitable due to the multiplicative loss of any infor- 
mation transferred recursively along the chain [3|. 

We demonstrate here that this assumption is in fact in- 
correct and that unexpected long range correlations re- 
main. They are responsible for a scale free power law 
regime with Pi{s) = Ca{p)s~^ for g{p) <^ s N {g{p) 
being the number of monomers per blob at monomer den- 
sity p) characterized by an exponent ui > 1 and a density 
dependent amplitude. Our simulation results are pre- 
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FIG. 1: Sketch of the bond-bond correlation {In ■ lm=n+s) of 
the bond vectors h = rt+i — ri, ri being the position vector of 
monomer i. The bond-bond correlation function is defined as 
the first Legendre polynomial Pi{s) = (/„ • 'n+s)„ /l^ over all 
possible bond pairs with curvilinear distance s — m — n > 
normalized by the mean squared bond length l^ = {in}- The 
typical size of the corresponding chain segments is R{s) = 

{i~nm.Y ^ with r„m = rm^n+a — rn- The dashed lines show 
the four relevant graphs of the analytical perturbation calcu- 
lation. The numerical factors indicate the relative weights of 
the leading l/s''^^ contributions to Pi{s). 



sented first and discussed together with simple scaling ar- 
guments. We focus on the scale free limit and finite chain 
size effects (s — > N) are considered more briefly. The an- 
alytical calculation (summarized graphically in Fig.^ is 
outlined at the end of paper. 

Two standard simulation methods for coarse-grained 
polymer chains 7] have been used to equilibrate long 
flexible polymer chains for which the bond-bond corre- 
lation functions presented in Figs. [3 El and 01 below have 
been computed. The body of our data comes from the 
"bond fluctuation model" (BFM) — a lattice Monte Carlo 
scheme where a monomer occupies 8 lattice sites (i.e., the 
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volume fraction is 8p) and the bonds U between adjacent 
monomers can vary in length and direction, subject only 
to excluded volume constraints and entanglement restric- 
tions Using a mixture of local, slithering snake and 
double-bridging moves 0, have created ensembles 
with chain lengths ranging up to TV = 2048 for densities 
between p = 0.00195/8 to p — 0.5/8 contained in large 
periodic boxes of linear size L = 512 ^ Re which allows 
us to eliminate finite-size effects. As can be seen in Fig.|21 
we have also studied single chains up to iV = 32768 in an 
infinite non-periodic box (p = 0) to characterize properly 
the dilute reference point. Additional molecular dynam- 
ics simulations of the bead-spring model discussed in 
have been performed to check explicitly that our results 
are not caused by lattice artifacts. As shown in Fig. 01 
for one example with N = 256 and p = 0.83, identi- 
cal behavior (although with reduced statistics) has been 
found. 

Bond-bond correlation functions for different densities 
are presented in Fig. [51 where we focus first on the be- 
havior of long chains (s ^ N). In the dilute limit 
{p 0) the athermal BFM chains are well fitted by 
R{s) = bos'^, V = 0.588 being Flory's good solvent 
exponent and &o ~ 3.0. The power law slope with 
w = Wo = 2 - 2z/ « 0.824 and Ca = v{2v - l){bollf 
observed in this limit (dashed line) is in fact expected 
[lo| from the formula 

{In ■ Ini) = (9„r„ • dmr,n) = -^9„9m {rnyn) (1) 

expressing Pi{s) as the second derivative of R{s). It is 
crucial that 21/ > 1 and, hence, < Ij i-e. the in- 
tegral over the correlation function is dominated by its 
upper bound. Otherwise, lo could not be governed by 
the exponent v describing the asymptotic chain size. In- 
terestingly, it follows also directly from Eq. for ideal 
gaussian chains (2j/ — > 1 : Ca — > 0) without additive power 
law corrections that Pi (s) must decrease stronger than 
any power law, i.e. exponentially. 

Coming back to Fig. (21 we note that also for finite den- 
sities Pi{s) coincides with the dilute power law for small 
curvilinear distances s where each chain segment inter- 
acts primarily with itself. At larger s ^ g{p) where the 
chains overlap and form a "melt of blobs" the correla- 
tion function decreases much faster, however not expo- 
nentially as one might expect, but with a second power 
law regime with w « 3/2 (bold lines). For the BFM the 
density p = 0.5/8 corresponds to a polymer melt since 
the excluded volume size becomes of order of the average 
bond size 0. We observe indeed that the dilute slope 
disappears in the melt limit and the correlation function 
can be fitted (bottom bold line) over nearly two orders 
of magnitude by 



Pl{s) = CaS-^ with UJ = 3/2, Ca = ^ ' (2) 
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FIG. 2: The bond-bond correlation function Pi(s) for BFM 
systems of different densities as indicated in the legend. The 
two lines indicate power laws corresponding to the asymp- 
totic regime for dilute {loq « 0.824, dashed line) and dense 
(uj — 3/2, bold lines) solutions, respectively. The observation 
of the second power law regime is the central result of this 
work. For systems in the so-called semidilute regime where 
the excluded volume is sufficiently large both exponents can 
be seen. The curvilinear distance at the crossover corresponds 
to the number of monomers per blob g{p). 

where we anticipate the analytical prediction for the 
asymptotic scale free limit. Considering that all param- 
eters are known the agreement is excellent although a 
slightly better fit could be obtained by accepting a larger 
exponent. Obviously, the chain statistics in the semidi- 
lute and melt regime must become in leading order gaus- 
sian, i.e. the integral over Eq. must be dominated by 
the lower bound in agreement with the finding uj > 1. In 
fact, Eq. is consistent with 

Risf^bh- s'- (3) 

[LU — 1)(2 — UJ) 

with b being the statistical segment length As long 

as LO < 2, the correction term (being due to an upper 
integration bound) contributes a negative, but ever de- 
creasing contribution to the chain distance, i.e., i?,(s)^/s 
approaches b^ from below. Interestingly, Eq. Q shows 
that the computation of Pi{s) allows to highlight directly 
corrections to gaussian behavior which can be blurred by 
the large ideal term if moments of spatial distances are 
considered. 

The scaling of Pi{s) and with density are analyzed 
in Fig.O Following the classical density crossover scaling 
U we take as natural scaling variable x = s/g{p) and 

Ms)9ipr - m - {::: (4) 

This is well borne out by the data. To be specific, we 
have used g{p) oc although the power law 

form cannot hold strictly for large densities where the 
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FIG. 3: Crossover scaling for Pi(s) for the BFM (TV = 2048) 
for different densities. In the main panel, the rescaled bond- 
bond correlation function f{x) = Pi{s)g{p)^° is plotted ver- 
sus the natural scaling variable x = s/g{p) assuming g{p) oc 
^-i/(3i— 1) j^Q^g ^^^^ £jjj. jj^j.gg g^pj £j^g gjjg^j cut-off (see 

Eq. becomes visible. In the inset we explicitly verify 
the density dependence of the amplitude Ca (p) obtained from 
Pi(s) (spheres) and R{s) (squares) using Eq. @ and Eq. 
respectively. The prediction for the melt regime, Eq. is 
indicated by the dashed line, the scaling Ca ~ ff(p)'^~"'' ~ 
^-0.885 £qj. £jjg semidilute regime by the bold line. 

semidilute blobs become too small and the assumed den- 
sity dependence must be a crude estimate. Considering 
this, the scaling collapse for all densities (obtained with- 
out any arbitrary shift parameter) is far from obvious 
and quite satisfactory. The amplitude Ca{p) of the non- 
gaussian corrections has been fitted directly from Pi{s) 
and R{s) assuming lu — 3/2. As can be seen in the inset, 
the data compare rather well with the prediction for the 
semidilute regime Cq ^ g{p)'^^'^° ^ p^O-^ss^ 

In the following we concentrate on dense melts. In 
Fig. ^ chain length effects are discussed at constant den- 
sity. Data from the bead-spring model have also been in- 
cluded to demonstrate the universality of the result. To 
collapse this correlation function on the BFM data it has 
been vertically shifted by the ratio of the amplitudes Ca 
calculated from Eq. ^ for BFM and bead-spring model. 
As can be seen for iV 16, exponentials are compatible 
with the data of short chains. (This might explain how 
the power law scaling has been overlooked in previous 
numerical studies since good statistics for large chains 
{N > 1000) has only become available recently.) How- 
ever, it is clearly shown that Pi{s) approaches systemat- 
ically the scale free asymptote with increasing N. The 
departure from this limit is fully accounted for by the the- 
ory if chain end effects are carefully considered (dashed 
lines). Generalizing Eq. perturbation theory yields 

Pi (s) ^ CaS-^/^a{x) + CbS-^/^b{x) (5) 

where we have set x — ^ s/N and Cf, — Aca/N and 




FIG. 4: Bond-bond correlation function versus s at melt den- 
sities for different iV for the BFM (p = 0.5/8) and the bead- 
spring model (p — 0.83) for A*' — 256. The bold line indicates 
Eq. ^ with Ca « 0.09 for long BFM chains. The dashed- 
dotted curve Pi(s) ~ exp( — s/1.5) shows that exponential be- 
havior is compatible with the small A*' data. The dashed lines 
correspond to the complete theoretical prediction Eq. Q for 
finite chain lengths A'' = 16, 32, 64 and 256 (from left to right). 
Considering that the theory does not allow for any free fitting 
parameter the agreement is good. In the inset we check di- 
rectly the recursion relation Eq. (|SJ by plotting K{s) versus 
U{s) — s/R{s)^p (same symbols as in main figure). The pro- 
portionality of K{s) and U{s) is well confirmed (bond line) 
by our data for 1 <^ s <^ N . 



defined the functions a{x) — b{x){l + 3a; + x'^) and 
b(x) = (1 - x)^/{l + x). For a; < 1 only the first 
term contributes: however, as s — > both terms be- 
come of similar magnitude. Both contributions to the 
correlation function vanish rigorously in this limit, and 
Pi(s)«(l-x)2. 

We outline now very briefly how Eqs. JSJ and 
have been obtained using perturbation calculation |l3| . 
Closely following Edwards [3 we first determine (r^^) « 

{rlm)o (1 + (^)o) - {TlmU)o and then, using Eq. 
the bond-bond correlations. Here, (...)o denotes the av- 
erage for the distribution function of the unperturbed 
ideal chain and U = dk dl Vi {rj^i ) the effective 
perturbation potential. The simplest approximation 
for the pair interaction potential in Fourier space is 
■^1(9) = ^(1^ l{<f + ■with V being the bare excluded 
volume parameter and ^ = (Z^/12pu)^/^ the screening 
length. This is sufficient for calculating the scale free 
regime (Eq. Q). The graphs which contribute to {In ■ Im) 
and their relative weights are indicated in Fig. ^ Note 
that the interactions described by the strongest graph Iq 
align the bonds In and Im while the others tend to reduce 
the effect. More care is needed to describe properly the 
finite chain size corrections Eq. (O , and the Pade approx- 
imation for the intrachain structure factor must be used 
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here. It can be readily shown that the different integrals 
are always dominated by long wavelength physics where 
Vi{q) — > q^P/12p ~ v'^. Therefore, the coefficients and 
Cf, of Eqs. (|2Jl and Q do not depend on local interaction 
parameters, such like v Q|. 

Since these calculations are lengthy we present a simple 
scaling argument for the uj exponent which captures the 
central physical idea (31] ■ As suggested by Eq. Q a direct 
measure for the non-gaussian corrections can be defined 
hy K = (i?2(2s) - 2i?2(s))/(2i?2(s)) ~ CaS^^'^ where 
we compare the size of a segment of length 2s with the 
size of two segments of length s joint together. Equiv- 
alently, this can be read as a measure for the swelling 
of a chain where initially the interaction energy U be- 
tween the two halves has been switched off. Obviously, 
K must be a functional of U and K{U = 0) = 0. Per- 
turbation theory tells us that the lowest order of the ex- 
pansion of K in U must be linear. What is still missing 
is an estimate for the scaling of the interaction energy 
between the segments. For incompressible melts it is well 
known that the energy penalty for joining two chains 
of length N is set by the depth of the correlation hole, 
U{N) « N/R{NY/p ~ N^-'^" This generalizes 

to arbitrary s. Hence, in c? = 3 the interaction energy is 
small {U{s) ~ which finally justifies the pertur- 

bation calculation. Taking everything together we obtain 
Lo = dv = ^1/2 and the important recursion relation 

The recursion relation provides a very compact and en- 
tirely self-consistent description of dense polymer melts 
(d > 2). To lowest order it leads to Eqs. I^Jl and ||2Jl. This 
shows that while the chain configurations must be gaus- 
sian to leading order, corrections do necessarily occur — 
even for very long chains. Eq. (jSJ has been directly vali- 
dated in the inset of Fig. ^ where we have plotted K{s) 
versus U (s) . The predicted linearity is well confirmed 
for large segments (s ^ 1). The recursion relation be- 
ing independent of the total chain length N does not 
capture the scaling for very large segments. Hence, the 
breakdown of the scale free behavior for K{U) (inset) and 
Pi (s) (main figure) are both due to the same finite chain 
size effect. 

In summary, we have shown that long range correla- 
tions exist in polymer melts and semidilute solutions due 
to both chain connectivity and incompressibility. Their 
most striking effect is the power law asymptote for the 
bond-bond correlation function allowing a direct numer- 
ical test of the non-gaussian corrections suggested by 
the self-consistent recursion relation Eq. ©. An im- 
portant consequence of this work arises for an exper- 
imentally relevant quantity, the static structure factor 
S{q). In fact, simulation and theory show distinct non- 



monotoneous behavior of q^S{q) vs. q (Kratky plot) due 
to the non-gaussian corrections. This suggests a possible 
route for experimental verification and is cause of serious 
concerns with regard to the standard operational defi- 
nition and measure of the persistence length from the 
assumed Kratky plateau. Finally, we point out that the 
physical mechanism which has been sketched above is 
rather general and should not be altered by details such 
as a finite persistence length — at least not as long as 
nematic ordering remains negligible. This is in fact con- 
firmed by preliminary and on-going simulations. 
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